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Farming and herding were introduced to Europe from the Near East and Anatolia; there are, however, considerable 
arguments about the mechanisms of this transition. Were it the people who moved and either outplaced or admixed 
with the indigenous hunter-gatherer groups? Or was it material and information that moved — the Neolithic Package — 
consisting of domesticated plants and animals and the knowledge of their use? The latter process is commonly referred to 
as cultural diffusion and the former as demic diffusion. Despite continuous and partly combined efforts by archaeologists, 
anthropologists, linguists, palaeontologists and geneticists, a final resolution of the debate has not yet been reached. In 
the present contribution we interpret results from the Global Land Use and technological Evolution Simulator (GLUES). 
This mathematical model simulates regional sociocultural development embedded in the geoenvironmental context during 
the Holocene. We demonstrate that the model is able to realistically hindcast the expansion speed and the inhomogeneous 
space-time evolution of the transition to agropastoralism in western Eurasia. In contrast to models that do not resolve 
endogenous sociocultural dynamics, our model describes and explains how and why the Neolithic advanced in stages. We 
uncouple the mechanisms of migration and information exchange and also of migration and the spread of agropastoralism. 
We find that (1) an indigenous form of agropastoralism could well have arisen in certain Mediterranean landscapes but 
not in northern and central Europe, where it depended on imported technology and material; (2) both demic diffusion by 
migration and cultural diffusion by trade may explain the western European transition equally well; (3) migrating farmers 
apparently contribute less than local adopters to the establishment of agropastoralism. Our study thus underlines the 
importance of adoption of introduced technologies and economies by resident foragers. 
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1. Introduction 

The transition to agropastoralism in western Eurasia 
between 10 000 and 3000calBC was associated with 
enormous cultural, technological and sociopolitical 
changes. Growing crops and herding animals have 
profoundly changed and continue to change global human 
history (e.g. |Roth| 



1.1. Archaeology 

Ncolithisation is believed to have begun during the early 
Holocene in the Fertile Crescent, a mountainous region 



between the Levantine coast and the Zagros ridge (Flan 



nery 



19731. Archaeobotanical, archaeozoological and ar- 



et al. 



2010 



1887 



Mithen 20041 Barker 2006 Ruddiman 



Westropp 



1872 



Diamond 



2002 



2006 Kaplan 



Kutzbach et al. 2010). These changes may 



be viewed positively as a trajectory of progress in the 
way it had been seen by nineteenth and early twentieth 



chaeogenetic work has demonstrated that all food crops 
and animals — except the dog — have their origins in and 
around the Fertile Crescent as a single founder region. 
The assemblage making up the Neolithic Package includes 
wheat, barley, rye, lentils, peas ( Willcox| 2005), and cattle, 



century evolutionists (Westropp 1872 Childe 19361, or it 



2007 Larson et al. , 2007 Zeder 2008 ) 



sheep, goat and pigs (Luikart et al. 2001 Edwards et al. 



may be seen as a road to perdition as it was, for instance, 
considered by J. Diamond (1997): "... a catastrophe 
from which we have never recovered. With agriculture 
came the gross social and sexual inequality." However 
the interpretation, the transition to agropastoralism, 
often termed Neolithisation, constitutes a major period of 
change in the history of humankind. 



While tendencies towards sedentism and storage of wild 
plants may already be interpreted from the Natufian 



data (Boyd 2006), intensive cultivation and domestica- 



tion of both plants and animals gradually began during the 
Younger Dryas and only fully developed during the early 
Holocene ( |Zeder| |2008| |Willcox et "aL| |2009| ). European 
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agropastoralism is allochthonous and its most likely ori- 
gins are in the Fertile Crescent, where farming and herding 
began — still in mixture with a broad spectrum of foraging 
practices — during the tenth millennium cal BC ( Flannery 



1973 Kuijt and Goring-Morris 2002). 
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The wider expansion of agropastoralism started around 
8500calBC, approximately 1000 years after the first ap- 
pearance of domesticated cereals in the Levant. The first 
clear evidence for colonist farmers was found on Cyprus 



In a later publication, IChildel (119421) suggested that popu 



(Peltenburg et al. 


2000 


Colledge et al. 




2004 


Willcox 


2005 


1; the expansion ended after 4000 ca 


BC, when Neo- 



lithic sites emerged on the British isles and throughout 



northern Europe (Sheridan 2007 Whittle 2007). Details 



of the intermediate region specific accounts of transitions 



have been collected, for example, by Price (2000), Whittle 



(2007), and Gronenborn and Petrasch (2010), including 
the prominent sixth millennium linear pottery cultures of 
central Europe (LBK, e.g. Liming and Stehli 19941 and 



the funnel beaker culture of the northern European plains 



(TRB, after 4500 calBC, |Midgeley[ |1992| . Not only did 
agropastoralism spread to the northwest from the Near 
East centres but also eastward as far as the Indus valley 



(Fuller 2006 



The question as to why agropastoral life style spread has 
been recently connected to environmental variations and 
conflict resolution 



Dolukhanov (1973), Weninger et al 
(2009]) and Gronenborn (2009a I suggested an emergence 



and spread of farming as a result of climate induced crises 
periods, during which it may have become necessary for 
groups to fission, i.e. to move from one location to the 
other to escape conflicts. 

Two contrasting concepts on the mechanism of the 
spread of farming across western Eurasia have existed 
side by side. One suggests the introduction of the new 
agropastoral technologies through movements of people — 
migrations of any form; the other suggests a techno- 
logy shift through indigenous adaptations and inven- 
tions fostered by culture contact — information dispersal 



of any form. Zvelebil ( 1998 ) discriminates seven spreading 
modes, for example elite exchange or leap-frog colonisa- 
tion, as combinations or intermediate forms of the two 
opposite spread mechanisms. 

The acculturation or cultural diffusion model corpus 
has, in the more recent past, been applied by a num- 
ber of post-processual archaeologists — more typically for 



the British isles but also for the continent (Hodder 1990 



Thomas 1991). It may go back to a critique by, for ex- 



ample, Zvelebil and Zvelebil ( 1988 ) of a migrationist model 



proposed by Renfrew ( 1987) and later by others (e.g. Bell- 
wood 2005 ). In a way connected to these models are those 
where farming or animal husbandry were seen as regional 
developments within Europe. Such indigenist scenarios 
have been proposed for southern France, northwest Africa 



or Greece (e.g. 


Geddes 


1980 


Courtin and Erroux 


Winiger 


1998 Theocharis and Bokonyi 


1973 


)■ 



Opposing this position is the one of migration, where 
the new technology and cultigens arrive from Anatolia 
and the Fertile Crescent into Europe through migrating 
people. This migrationist position goes back to [1925 , 
when V. G. Childe| noticed a gradient in the spatiotem- 
poral distribution of ceramics from western Eurasia em- 
anating from the Fertile Crescent northwest into Europe. 



lation pressure in the source region was the driver of this 
outmigration. His position was supported by |Ammerman 
and Cavalli-Sforza ( 1971 1973 ), who formulated the 'wave 



of advance model', which was based on the concept of de- 
mic diffusion. Today, scholars from a number of disciplines 
favour migrationist models both for people and for culti- 



gens (Sokal et al. 


1991 Richards| 2003 


Pinhasi et al. 


2005 |Edwards et a 


l.[ 2007| Bramanti et a 


l.| |2009| Balar- 


esque et al. 


2010 


Haak et al. 


2010 


1. Archaeology, partic- 



ularly continental European archaeology, sees evidence for 
more complex scenarios of migrations and local accultur- 
ation notably in western central Europe and France (e.g. 
Jeunessel 120001 iBentley et al.l 120021 iGronenbornl |2007bb, 



but acknowledges that long-distance contacts across west- 
ern Eurasia did exist during the mid-Holocene and should, 
at least partly, have been maintained by the migration of 
people. It is yet unclear exactly when migrations began 
and what the relative importance of acculturation and 
movement of people was. It may not be ruled out that 
the large-scale population replacements around the Neo- 
lithic began with the onset of the sixth millennium cal BC 
( Gronenbornl |2007al [20TT ) . 



1.2. Mathematical Models 

The spatiotemporal structure of the advance of farm- 
ing in Europe was first and very coarsely quantified by 
Edmonson (1961), who estimated the speed of the agro 
pastoral transmission frontier at 1.9 km a -1 . Later, 
(119651) and lAmmerman and Cavalli-Sforza ( 1971 



Clark 



19731 



based their analysis on radiocarbon dates (in three areas, 
at 53 sites, 103 sites, respectively) to calculate the velocity 
of the appearance of farming practice in Europe along a 
southeast-northwest gradient. All three studies found an 
approximately linear relationship between temporal and 
spatial distance of European Neolithic sites to four Near 
Eastern sites, with a slope of approximately 1 km a - 1 . Pin- 



hasi et al. (2005) confirmed this finding on a more exten- 
ded data set of 765 sites; they calculated a spreading rate 
between 0.8 and 1.3 km a -1 ; even when shortest-path dis- 
tances are considered (longer than great circle distances 
because of the detour necessary from the Near East to 
Anatolia and from central Europe to Iberia), a similar rate 
(0.6-1.1 kma -1 ) is found. A data set of 477 sites, including 



boreal European sites, was used bylDavison et al. (2006 



2007| 20091 who simulate for a nonintcracting agropas- 
toral subsistence style Neolithisation and its speed with a 
reaction-diffusion model. They too arrive at a mean speed 



1974 of lkm a into Europe. 



Ackland et al. (20071 simulated the spread of farming 
by including a 'hitchhiking' advantageous trait in their 
reaction-diffusion model. This new lifestyle addition could 
either reflect the immigration of farmers, or a dynamical 
conversion of foragers into farmers upon contact. Their 
model predicts that Neolithic farmers outplaced the indi- 
genous population up to a line approximately connecting 
today's Venice-Prague-Warsaw-Moscow. To the north 
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Figure 1: Map of Europe with simulation region boundaries (solid 
lines) in the Global Land Use and technological Evolution Simulator 
(GLUES). Dots represent 631 radiocarbon dated Neolithic sites from 
|Pinhasi et al,| ( |2005| ) ■ Our discussion focusses on a southeast to 
northwest transect along the highlighted regions A to J. 



and west of this line, the so-called converts adopted the 
new lifestyle by cultural diffusion. 

Available evidence, in part originating from isotopic and 
genetic studies, points to a discontinuous expansion se- 



quence for western Eurasia ( 


Guilaine 


2001 


Gronenborn 


2009a 


Bocquet-Appel et al. 2009 Schier 


20091 during 



which short dynamic phases of long distance rapid expan- 
sions were followed by periods of stand-still with local or 
regional colonisation. Discontinuities in the Neolithic ad- 
vance, however, have not been hindcasted by the aforemen- 
tioned models. One possible reason is that in the frame- 

(e-g- 



works provided by Ackland et al 



et al. 



20071 or Davison 



These 



(2006 20091 important aspects are missing 
may comprise more detailed descriptions of the resources 
needed and used by the people, the influence of the local 
biogeographic suitability for farming or herding, and of 
temporal variation in resource availability. Their mod- 
els do not simulate any endogenous cultural, technical or 
agrarian development. All these factors may in principle 
accelerate or slow down the process of Neolithisation and 
lead to a more complex spatiotemporal pattern than may 
be predicted by simple reaction-diffusion models. 

In this study, we employ the Global Land Use and tech- 
nological Evolution Simulator (GLUES), which resolves 
local innovation, migration and cultural diffusion of traits 

Lemmen 2009| Lemmen and 



(Wirtz and Lemmen 2003 



Wirtz| [2010J). Although GLUES has been developed for 
the global domain, we restrict our analysis in this study to 
western Eurasia, where radiocarbon dates from Neolithic 
sites are abundant and of high quality, and where the is- 
sue of migration versus cultural diffusion is most intens- 
ively debated. The model is chosen because it allows us to 
differentiate between exchange (i.e. information exchange 
as cultural diffusion) and migration (demic diffusion) as 
important vectors of the expansion of agriculture. 



In the following section, we shortly introduce the 
GLUES model and the radiocarbon site data which are 
used for validation; a full description of the algorithms 
used in GLUES can be found in the supplementary on- 
line material (SOM). The spatiotemporal pattern of the 
emergence and advancement of agropastoralism in west- 
ern Eurasia is reconstructed and analysed in detail. This 
is achieved through a model-data comparison for ten focus 
regions along a southeast-northwest trajectory, from the 
Levant to north Germany; model-based expansion rates 
are put into the context of prior estimates from radiocar- 
bon dates. A major part of our discussion concentrates on 
the discrimination of migration versus trade and the max- 
imum contribution of immigrants to emerging agropastoral 
communities in Europe. 

2. Material and Methods 

GLUES mathematically resolves the dynamics of local 
human populations' density and characteristic sociocul- 
tural traits in the context of a changing biogeographical en- 
vironment. A local sociocultural coevolution is described 
by changes in mean population density, technology, share 
of agropastoral activities, and economic diversity, within 
a simulation region of approximately country-size extent 
(Figure [I]). Each local population utilises its regional 
natural resources, which are described by vegetation pro- 
ductivity and climatic constraints. Each local population 
interacts with its geographical neighbours via trade and 
migration. The conceptual model is outlined below, for 
details on the algorithms used and the mathematical im- 



plementation we refer the reader to Wirtz and Lemmen 



(2003), summarised in the SOM. 



2.1. Characteristic traits 

For pre-industrial human societies, we define three char- 
acteristic traits: 

1. Technology is a trait which describes the efficiency 
of food procurement — related to both foraging and 
farming — and improvements in health care. In partic- 
ular, technology as a model describes the availability 
of tools, weapons, and transport or storage facilities. 
It aggregates over various relevant characteristics of 
early societies and also represents social aspects re- 
lated to work organisation and knowledge manage- 
ment. It quantifies improved efficiency of subsistence, 
which is often connected to social and technological 
modifications that run in parallel. An example is the 
technical and societal skill of writing as a means for 
cultural storage and administration, with the latter 
acting as a organisational lubricant for food procure- 
ment and its optimal allocation in space and among 
social groups. 

2. A second model variable represents the share of farm- 
ing and herding activities, encompassing both animal 
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husbandry and plant cultivation. It describes the al- 
location of energy, time, or manpower to agropastoral- 
ism with respect to the total food sector. We define a 
local population as Neolithic when this share is larger 
than the share of foragers — regardless of its techno- 
logy, economic diversity, or population density. 

3. Economic diversity resolves the number of different 
agropastoral economies available to a regional popu- 
lation. This trait is closely tied to regional vegeta- 
tion resources and climate constraints. We do not, 
however, attribute specific plants and animals to each 
economy. As an example, a value of four would be ob- 
tained when (1) domestic pigs and (2) goats and the 
growing of (3) barley and (4) wheat were present in a 
given population. A larger economic diversity offering 
different niches for agricultural or pastoral practices 
enhances the reliability of subsistence and the efficacy 
in exploiting heterogeneous landscapes. 

The temporal change of each of these characteristic 
traits follows the direction of increased benefit for suc- 
cess (i.e. growth) of its associated population; this concept 
had been derived for genetic traits in the works of Fisher 
(1930), and was recently more stringently formulated by 



Metz and colleagues (Metz et al. 1992 Dieckmann and 



Law 1996 Kisdi 2010) as adaptive dynamics (AD). In 
AD, the population averaged value of a trait changes at 
a rate which is proportional to the gradient of the fitness 
function evaluated at the mean trait value. The AD ap- 
proach was extended to functional traits of ecological com- 



munities ( 


Wirtz and Eckhardt 


1996 


Merico et al. 


2009 


Smith et al. 


2011 


), and was first applied to cultural traits 



Lcmmcn (2003). 



of human communities by Lemmcn (2001) and Wirtz and 



2.2. Local 



resources 



Each simulation region is defined by a largely homo- 
geneous vegetation productivity (measured as net primary 
productivity, NPP), resulting in an average size of 130 • 
10 3 km 2 (Figure [TJ. We reconstruct past distributions of 
NPP with a global climate model coupled to a vegeta- 



tion module. Climber-2 ( Claussen et al. 1999 1 tempcrat 



ure and precipitation anomalies from the IIASA climato- 
logical data base (International Institute for Applied Sys- 



tems Analysis, Leemans and Cramer 1991 1 are converted 



to NPP according to the climate constraints on NPP from 



Lieth (1975); we do not use soil maps to constrain veget- 
ation productivity. 

From NPP, both the regional utility of natural food re- 
sources and the number of potential domesticates are de- 
rived. According to |Braidwood and Braidwood s ( 1949 



1950) hilly flanks hypothesis, potential domesticates were 
most abundant in open woodlands at low to intermediate 
NPP. The number of potential domesticates furthermore 
depends on a continental aggregation to account for the 



2.3. Exchange of information and people between regions 

Information exchange and migration are vectors of the 
spread of technology, economic diversity, and farming 
practice from the founding centres to adjacent simulation 
regions. We discriminate the diffusion of traits without 
involving resettlement of people (cultural diffusion by in- 
formation exchange), and the diffusion of traits via mi- 
gration (demic diffusion). In GLUES, both mechanisms 
are driven by differences in influence between neighbour- 
ing local populations. 

We assume that information travels two orders of mag- 
nitude faster than people. Exchange networks extend 
over distances of up to 1000 km, in the later Mesolithic 



and Neolithic (Mauvilly et al. 2008 Gronenborn 1999); 



these networks were crossed many times during the active 
time — say ten years — of a Neolithic trader. Within this 
time span, a migration model like the one by | Ammerman 
and Cavalli-Sforza] (19731, would allow for an advance of 
only ten km. This parameterisation leads to diffusivities 
for migration on the order of ten km Si . Si value which 
is comparable to the diffusivities employed by other model 



studies of demic diffusion (e.g. Davison et al. 


2000 


Ack- 


land et al. 


2007 


Patterson et al. 


20101. 



2.4- Reference data and simulated time scale 

Our reference data set is the comprehensive data col- 
lection of 765 sites by Pinhasi et al. ( 2005 1 . These au- 



area-biodiversity relationship (e.g. Begon et al. 19931. 



thors used site data provided by the United Kingdom Ar- 
chaeology Data Service, the Central Anatolian Neolithic 
e- Workshop (CANeW), the radiocarbon CONTEXT data- 
base, and the Radiokarbondaten Online (RADON) data- 
base. In their compilation, they included only sites with 
small dating uncertainty (< 200 a); they report dates as 
calibrated calendar years before present (relative to 1950) 
based on calibration of original 14 C measurements with 
CalPal 2004. This data set was created by the |Pin hasi 
|et al.| to provide a high quantity of dates and good spa- 
tial coverage at the expense of chronologic uncertainties, 
which could have been avoided, if for example only AMS- 
dated (accelerator mass spectroscopy) samples had been 
used; only few of the 765 sites, however, have been AMS 
dated. For our purpose, this data set with many (possibly 
uncertain) dates represents the expansion of agropastoral- 
ism at a satisfactory level of detail. Future simulations at 
a refined spatial scale would benefit from a data set with 
better chronologic control, where local and regional events 
are presented in higher resolution and where the region- 
ally patchy nature of the expansion of agropastoralism is 
better represented. 

From this data set, we choose for comparison those 
631 sites which are located in the spatial model domain 
(10°W-42°E and 31°N-57°N) and the period of interest 
(8000-3500 calBC). For each site, we use the age range 
computed from the reported calibrated radiocarbon age 
and the reported standard deviation. 
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For the mathematical model, we introduce the age scale 
'simulated time BO (simBC) to distinguish between em- 
pirically determined age models and the model time scale. 
Ideally, simBC should be numerically equal to calBC. 

We set up the eight global model parameters such that 
the simulation is able to hindcast an accurate timing 
and location of the early farming centres Fertile Crescent, 
northern China, and Mesoamerica (Smith 1997), and a 



reasonable global pattern of the subsequent Neolithisation. 
The simulation is started at 9500 simBC. All of the 685 
biogeographically defined regions (including 71 in western 
Eurasia) are initially set with farming activity at 4% and 
established agropastoral communities at 0.25, what repres- 
ents a low density Mesolithic technology population and 
a broad spectrum foraging lifestyle with low unintentional 
farming activity. The latter is assumed to represent early 
animal harvesting, selective seed gathering, and the active 
use of fire. 



3. Results 

In the GLUES simulation, farming originates in the Le- 
vant (focus region A, cmp. Figure ^ around 7000 sim BC 
and penetrates into Europe in a northwest direction. By 
3500 simBC all of continental Europe has converted to 
farming as the predominant subsistence style. This emer- 
gence of farming in western Asia and Europe is shown as a 
series of snapshots in Figure [2] (see SOM for an animated 
version with finer temporal resolution). 

3. 1 . Expansion of agropastoralism 

The initial development progresses slowly and at a low 
level. It begins during the first century of the seventh 
millennium simBC in a region encompassing today's Le- 
banon, coastal Syria and a small part of the adjacent 
coastal Anatolia. In the 67th century simBC, northern 
Greece converts to agropastoral subsistence with rapid 
extension into the central Balkan. Over the next four 
hundred years, these agropastoral nuclei spread out fur- 
ther, encompassing the whole of Greece and the south- 
ern Balkan, and the coast of Anatolia by the 63rd cen- 
tury sim BC. 

A rapid expansion of agropastoralism occurs between 
6200 and 6000 simBC, transforming the entire Balkan re- 
gion and Anatolia. By 5750 simBC, the new subsistence 
mode has reached the northwestern and the easternmost 
coasts of the Black Sea. In the 57th century simBC, in- 
dependent agropastoralism arises in north Africa in the 
region around the Strait of Gibraltar, and it emerges on 
the Italian peninsula. 

The 55th century sim BC sees a rapid expansion of farm- 
ing and herding into the area of the central LBK, and its 
spread into the south coast of the Iberian peninsula. By 
the 54th century, the LBK has expanded west- and east- 
wards and covers a vast stretch of land from southern Ger- 
many to the Ukraine; this central-eastern European area 
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Figure 2: The spread of agropastoralism in western Eurasia from 
7500 simBC (top left, then downwards) until 3500 simBC, hindcas- 
ted with the Global Land Use and technological Evolution Simulator. 
A finer temporal resolution (50 a time step) animation of this evolu- 
tion is available in the SOM. 



intensifies agropastoral activity without notable expansion 
until 5100 simBC. 

Around 5000 sim BC, forager societies on the north coast 
of the Black Sea, in north Africa and on the Iberian penin- 
sula have converted to predominantly agropastoralism. At 
4750 simBC, the Neolithic package reaches the Baltic Sea 
at the Oder river mouth; this coastal agropastoralism ex- 
pands eastward until 4500 simBC and resembles the rise 
of the eastern TRB culture. By this time, farming and 
herding have — in the model — reached the south coast of 
France and the north coast of Portugal. 

An agropastoral area resembling the western TRB ap- 
pears by 4400 simBC, also in southern Germany the new 
life style becomes dominant. The later half of the fifth 
millennium sees a slow expansion towards the northeast 
of Europe, and the gap closure in central and northern 
France. After 4000 simBC, agropastoralism reaches the 
British isles. 
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Figure 3: Timing of the transition to agropastoralism in Western Eurasia. The simulated transition (background pastel shading) is contrasted 
with the radiocarbon ages of Neolithic sites from |Pinhasi et al.| | |2005| solid colour triangles). The lower right inset image shows the transition 
for a scenario without migration or exchange, i.e, it shows the propensity of regions to endogenously develop agropastoralism. 



3.2. Timing of agropastoralism in model and data 

A summary description of the timing of agropastoralism 
between 7500 and 3500simBC is illustrated by Figure [3] 
shown alongside are the median radiocarbon dates of Neo- 
lithic sites within this period from the data compilation by 



Pinhasi et al. (20051. From this time-integrated perspect- 



ive, the simulated centres of agropastoralism in the Fertile 
Crescent, in northern Greece and at the Strait of Gibral- 
tar are evident, as well as the southeast to northwest tem- 
poral gradient of the Neolithic transition. The model-data 
comparison shows many good matches between radiocar- 
bon dates and simulated transition dates. We can clearly 
see, however, the spatial scale difference between simula- 
tion region and site data. The spatial distribution of ra- 
diocarbon dated sites has good coverage along the transect 
from the Levant to northwestern Europe discussed below, 
it provides few or no information on eastern Europe, on 
the Iberian peninsula, and in north Africa. 

To assess the quality of the simulated onset of agropas- 
toralism, we compare in Figure [4] the change in fractional 
agropastoralism to the radiocarbon site statistics for ten 
focus regions A-H (a transect from the the Levant A to 
north Germany H) and radiocarbon dated sites within the 
region, or within 200 km distance of the region centre for 
small regions. We also indicate by colour selected cultural 
attributions. 

Seventeen sites within or near region A are dated 
between 8000 and 5500 BC, of which the most frequent 
cultural attribution is Pre Pottery Neolithic (PPNX, 
9 sites) and Pottery Neolithic (6 sites). The most fre- 



quent century is the 68thcalBC (4 sites); the simulated 
change in agropastoral activity is greatest in the 70th cen- 
tury simBC. 

Ten sites are found in or near region B, most of which are 
assigned to the Pottery Neolithic (6 sites) . All sites date to 
before 6000calBC, with a maximum around 6400calBC; 
the largest simulated change to farming occurs in this 
region and in region C, around 6300simBC. Near or in 
C, 15 sites cover a wide temporal range from 6500 to 
4400calBC. The site statistic within or around region D is 
poor with only six sites, which date to 5800-5000 calBC. 
The timing of the largest simulated change is 6100simBC; 
this simulated transformation resembles the occurrence of 
the Koros culture. 

Like region D, most sites near region E are attributed 
to Koros; the second most frequent cultural complex in re- 
gion E is the LBK, which is also the dominant attribution 
at sites around regions F to I. In regions H to J, the site his- 
togram is bimodal with the latter peak assigned to funnel 
beaker sites. Many site dates near region E fall within the 
period 6000-4800 calBC, whereas the simulated change is 
greatest at 6100simBC; around region F, the simulated 
transition occurs around 5800simBC. The most frequent 
date is the 52nd century calBC, with a large range of 1500 
years. 

Radiocarbon dates for region G range from 6600 to 
3500calBC; a maximum occurs between the 53rd and 
47th century, which is concurrent with the largest sim- 
ulated change at 5200simBC. Seven LBK sites around re- 
gion H are dated to 5600-5000 calBC, coterminous with 
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Figure 4: Timing of the transition to farming in GLUES and cal- 
ibrated radiocarbon date statistics for the ten focus regions (map 
inset). In each panel A— I (ordered from southeast to northwest, 
bottom to top), the change in the fraction of agropastoralism from 
GLUES (grey shading) is contrasted to the number of radiocarbon 
dates compiled by |Pinhasi et al.| ( |2005[ l (bars, including la uncer- 
tainties of the age determination). Colour indicates selected cul- 
tural attributions within the site record to pre-pottery Neolithic (of 
any kind, PPNJf), Koros, Linear pottery culture (LBK), and fun- 
nel beaker culture (TRB). To account for the scale difference, the 
simulated transition was broadened in time by a convolution with a 
normal distribution with mean zero and standard deviation propor- 
tional to region area and diffusion length. 



the simulated shift at 5500simBC. For region I, with 
mostly LBK-attributed sites and radiocarbon dates (5200- 
4400 cal BC), the simulated subsistence change culminates 
around 4600simBC. The respective model transition for 
region J appears at 4400simBC. Here, the site histogram 
can be divided into two modes, where the first encompasses 
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Figure 5: Agricultural onset timing versus distance from an assumed 
archaeological centre near Beirut (Lebanon) in the GLUES simula- 
tion and in the data set by |Pinhasi et ah] ( |2005^ - The correlation 
of distance and timing in the simulated data for n = 66 regions is 
r 2 = .39 (large circles, significant at p = .01), the average slope is 
0.81 km a" 1 . For n = 631 sites, the correlation is r 2 = .61 (blue dots, 
significant at p = .01), the slope is 0.72 km a -1 . The solid lines show 
for each distance interval of 500 km the probability of the earliest 
and farthest occurrence (p = 0.05 percentile) of agropastoralism in 
the site dates (blue) and the simulation (red). 



radiocarbon dates between 5800 and 4600 cal BC (incl. 9 
LBK sites) and the latter dates from 4200 cal BC (incl. 10 
TRB sites). 

3.3. Time lag-distance relationships 

The average speed of the expansion of agropastoralism 
can be estimated from the time lag-distance relationship 
relative to an assumed founding centre of agropastoral- 
ism (e.g. Ammerman and Cavalli-Sforza 1973). Figure [5] 



shows this relationship for all regions and radiocarbon 
dated sites within the model domain. Here, the assumed 
agropastoral centre is near today's Beirut, which lies in 
the middle of focus region A. The (great circle) distances 
and time differences to this assumed centre extend over 
4000 km and 4000 a, respectively. 

Time lag and distance from the sites are highly cor- 
related (r 2 = .61) but also indicate a stair-case like 



0.72 km a 1 (cinp. 



Guilaine 




2001 


Gronenborn 




2009a 



Schier 2009). This means that in the data collection the 
spread of Neolithisation is slower in spatial proximity of 
the founder region than is predicted by a linear correla- 
tion; between 6000 and 4500 cal BC, however, the graph of 
the majority of sites lies above the regression line, which 
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indicates a more rapid wave of advance from the Balkan 
towards central Europe. Lag and distance for GLUES- 
simulated regions are also correlated to a marked degree 
(r 2 = .40) and are similarly scattered around the regres- 
sion line. Of the ten focus regions, regions A and B develop 
more slowly than expected from the regression and regions 
E-G develop agriculture faster than the linear regression. 
The average speed for the expansion of agropastoralism 
from the Levant into Europe calculated from the model is 
0.81 km a" 1 . 



4. Discussion 

The Global Land Use and technological Evolution Simu- 
lator is able to hindcast a realistic spatiotemporal pattern 
of the introduction of farming and herding into Europe 
between 8000 and 3500simBC. The simulated expansion 
speed of agropastoralism compares well to a large dataset 
of radiocarbon dated Neolithic sites; the inhomogeneous 
spatial distribution of Neolithisation is reproduced. 

4-1- spatiotemporal onset and expansion of agriculture 

The differences we observe between simulated timing 
and the radiocarbon age of sites within a simulation re- 
gion (Figures [3] and [1]) are less than 1000 a for almost all 
sites, for the majority of sites less than 500 a; only a hand- 
ful of sites show differences greater than 1000 a. These 
differences are similar in magnitude to those obtained by 



Davison et al. (2007) between their numerical model and 
radiocarbon dated sites in Europe. At this scale of model 
uncertainty, the radiocarbon dating uncertainty of indi- 
vidual sites (< 200 a) can be neglected: the mismatch 
between the onset definitions in the data (presence of a 
Neolithic site, Figure |4| and in the model (50 % agropas- 
toral activity) , as well as the spatial scale mismatch (local 
site data versus country-size simulation region, Figure [3]) 
introduce larger temporal differences. To overcome the 
spatial scale problem, |Zimmermann| ( 2004 1 argued for a 
landscape approach to archaeology, whereby a multitude 
of local sites are used to infer the archaeological context 
at the regional scale or larger. The landscape approach 
can only succeed, however, if many sites within a region 
are excavated, as was the case for the lignite mining area 
of the Aldenhovener Platte studied by |Liining and Stehh] 
(1994) and Zimmermann (2004 1. From a model perspect- 



Gkiasta et 


al. (2003 r 2 


and Cavalli-Sforza ( 1971 


et al. ( 


2005 


r 2 = 0.64, n = 



= 0.53, n = 510), Ammerman 
r 2 = 0.79, n = 103), and|Pinhasi 



empirical results can be attributed to the number of sites 
under consideration, to the location of the assumed found- 
ing centre, to site selection, or to the consideration of the 



shortest land route versus great circle distance (Pinhasi 



et al. 20051. For calibrated dates, Pinhasi et al. calculate 
a speed range of 0.6-1.3 km a -1 when these differences are 
taken into account. The validity of comparing the onset 
of agropastoralism between simulated regions (with large 
areal extent) and (local) radiocarbon dated sites, despite 
the different scales, is supported by our simulation in two 
ways: (1) the marked correlation obtained between lag and 
distance of first agropastoralism to an assumed founding 
centre; (2) a calculated speed of 0.81 that agrees closely 
with other published estimates. 

Our simulations do not take coastal expansions into ac- 
count, which would seem a major model deficiency at first 
glance. The independent Moroccan model centre, how- 
ever, acts for the Iberian peninsula similar as an explicit 
fast migra tion process (like leapfrogg ing, on the order of 
20 km a" ' 



Zilhao 



1993 



Zilhao 



2000) along the Mediter- 



ranean coast and islands. GLUES does not currently ac- 
count for rivers: we attribute the late transition of north- 
ern France in the simulation to a missing pathway from the 
Mediterranean coast through the Rhone valley. Indeed, 



Davison et al. ( 2006 ) found in their model a significant role 
of waterways in the Neolithisation of Europe, whereas our 
results only indirectly (through the definition of regions 
by homogeneous vegetation) include river basins; GLUES 
performs well despite the lack of explicit river pathways 
like the Rhine or Danube valleys for all regions of Europe 
except central France. 

In the simulation, as well as in the data, the expansion 
of farming occurs in stages with periods of rapid spread 
followed by periods of local intensifications. The rapid 
Neolithisation from Greece to the central Balkan in the 
67th century simBC is followed by a 400 a period of rel- 
ative stagnation. A very similar pattern is hindcasted for 
the LBK-like Neolithisation in the 55th and 54th century 
calBC, and for the relative stagnation before the onset of a 
TRB-like Neolithic further north. Several regions exhibit a 
slow conversion to agropastoralism: according to Bellwood| 



and Oxenham s (2008) classification of zones in Neolithic 
Europe, France (see discussion of rivers above) would rep- 
resent rather a friction than a spread zone. Another fric- 
tion zone, where the Neolithic is introduced gradually, ex- 
ists in the northern European lowlands, where empirical 



data supports the simulated late arrival of farming (Midge- 



ley 


1992 


Zvelebil 


2006 


Hartz et al. 


2007) 



ive, more studies on methodologies to scale up the site 
(or many sites) information to the landscape are highly 
desirable. The scale difference illuminates the resolution 
limits of our model: GLUES resolves societal dynamics in 
larger environmental contexts rather than the history at 
individual sites. 

We find a marked correlation between the timing of first 
agropastoralism and the distance from a founding centre in 

the model (r 2 = .40), and an average speed of agropastor- spective, Davison et al. (2007. 2009) suggested that to 
alism in western Eurasia of 0.81 km a -1 . Using radiocar- 
bon data, a marked or high correlation was also found by 



Our simulations predict a second and eastern expansion 
path around the Black Sea, which was archaeologically 
suggested by Kotova (2003 2009). From a model per- 



find a suitable simulation consistent with the radiocarbon 
dated site context, one needed an additional early wave 
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of advance emanating around 8200calBC from an East- 
ern European centre; agropastoral expansion would then 
proceed via the steppe corridor and be responsible for the 
eastern version of the Neolithic in Europe. This steppe cor- 
ridor is archaeologically visible from the spread of pottery 
across Eurasia and the expansion of farming from east- 



ern Anatolia into the Caucasus (Dolukhanov et al. 2005 



Gronenborn 2009b Kotova 2003 2009). In the reaction 



diffusion migration model by Ackland et al. (2007), a cir- 



cular expansion from a single Mesopotamian centre is sim- 
ulated; the geographical bottlenecks between the Mediter- 
ranean and the Black Sea, and between the Black and the 
Caspian sea act as secondary wave centres and thus the ar- 
rival of the farming wave in Europe "appears to come from 
two sources: north and south of the Black Sea" (Ackland 



et al. p. 8715). GLUES exhibits the same behaviour, with 
a southern and eastern path around the Black Sea based 
on a single Mesopotamian source area; for central Europe, 
however, our model suggests that the secondary centre is 
rather located in northern Greece and the central Balkan. 
Archaeologically, the land route bottleneck may have not 
been too important for the Mediterranean coast, where the 
many islands provided a fast sea route from the Levant to 
Cyprus, Greece, and as far as Portugal ( |Peltenburg et alTj 
|Theocharis and Bokonylj |1973 Arias} |1999| |Zilhao 



2000 



2000). 



A separate non-Eurasian independent centre of agropas- 
toralism is simulated by GLUES in the Maghreb. From 
there, agropastoralism enters Europe via the Strait of 



Gibraltar around 5500simBC (cmp. Manen et al. 2007) 
Archaeological records to show this are sparse; it is clear 
however, that the strait had been in use as a migra- 
tion path long since pre-Neolithic times, which is evid- 



ent in gene pool analyses (e.g. Currat et al. 2010). Not 
only people but also domesticates crossed the strait in 
prehistory. This was verified in a study by |Anderung 



|et al. ( 2005 ) , who found mitochondrial DNA of Bronze age 
(1800calBC) Iberian cattle, of which a significant number 
possessed African haplotypes. 

4- .2. Demic or cultural diffusion 

The relative contribution of demic versus diffusive pro- 
cesses can be calculated by following the streams of mi- 
gration and trade in the model. We find that exchange 
processes contribute much less to local Neolithisation than 
adoption does. In Figure [6j we show this for the fraction 
of immigrant agropastoralists: for some (mostly Mediter- 
ranean) regions, immigrants are unimportant; for most re- 
gions, immigrants constitute one fourth of the agropastoral 
community; for a few northern and alpine regions, immig- 
rants dominate. 

We assessed the model sensitivity to different config- 
urations for the speed of both information exchange and 
migration. We find that (1) the model does not show sens- 
itivity to cither of these two parameters for a wide range of 
values, and that (2) the substitution fraction of agropastor- 
alists with immigrants generally remained around 25% for 




Percentage of immigrant farmers and herders 



Figure 6: Simulated fraction of immigrants in the agropastoralist 
population of each region at the time when the transition is loc- 
ally 90% complete. Green colour indicates mainly local adoption by 
resident foragers, yellow to red a major contribution of immigrants. 



most regions; local invention and adoption of ideas dom- 
inate the Neolithisation, irrespective of whether people or 
information moved. 

Is demic diffusion a sufficient explanation for the 
European Neolithic? This has been suggested by many 



authors (|Ammerman and Cavalli-Sforza 


1971 1973; Sokal 


et al. 


1991||Richards 2003 ; Edwards et al 


. 2007 


Bramanti 


et al. 


2009 


Balaresque et al. 


2010 


Haak et al. 


2010 


). We 



can confirm this finding from a model simulation where 
cultural diffusion was deactivated. Demic processes alone 
can reproduce the timing and lag-distance relationship 
seen in the radiocarbon data. 

But is demic diffusion necessary for explaining the ra- 
diocarbon record? Pinhasi et al. (2005) answered this 



question positively and pointed out that, at present, no 
working model existed that could explain the European 
Neolithic without demic diffusion. Alike the demic 
diffusion-only experiment, we set up a simulation experi- 
ment where only information was allowed to diffuse, and 
where migration was inhibited: we could successfully re- 
produce the spatiotemporal emergence of the Neolithic 
in Europe with purely cultural diffusion processes. To 
the question of necessity of demic diffusion for the ex- 
planation of the radiocarbon record, our answer is no. 
This shows that the recently published evidence for ma- 
jor population transfers around the time of the Neolith- 



isation process (Haak et al. 2010) may have to be func 



tionally disconnected from the spread of agropastoralism: 
apparently, people did move at greater scales during the 
sixth and fifth millennium but these movements were not 
triggered by the spread of farming. It may entirely be pos- 
sible that early — yet hypothetical — migrations already oc- 
curred before and during the seventh millennium and were 
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Region F 
agropastoralism 
migr:adopt 

demic 54:46 
mixed 22:78 
cultural 0:100 



Region H 
agropastoralism 
migr: adopt 

85:15 
41:59 
0:100 



Region H 
technology 
migr : exch : adopt 

22:0:78 
6:13:79 
0:21:79 



Table 1: Contributions (in percent) to local agropastoralism and 
technology from three different sources (1) demic diffusion (labelled 
migr), (2) cultural diffusion (exch), and (3) local adoption and in- 
vention (adopt) for three different model configurations with demic 
only, mixed (our standard configuration discussed in Figures 2-6), 
and cultural only diffusion. Simulation results are shown for model 
regions corresponding to today's Hungary (focus region F) and south- 
ern Poland (focus region H). 



undertaken by hunter-gatherers or mixed hunter-gatherer- 
horticulturalists originating from Anatolia. These immig- 
rants then gradually pushed the original Mesolithic hunter- 
gatherer population of Europe towards the continental 
margins. Later, this migratory stream was complemen- 
ted with farmers from Anatolia who then interacted with 



those hunter-gatherers who had arrived earlier (Gronen 



born 2011). This scenario would explain the archaeo- 



genetic evidence for migrations as well as the archaeolo- 
gical evidence of interactions by disconnecting the spread 
of farming from the mid-Holocene migratory processes. 

The insensitivity of the simulation results to the absence 
of either trade or migration processes prevents us from 
constraining the parameters for these processes quantitat- 
ively. Even more, it tells us that from the phenology (the 
timing of agropastoralism) we cannot infer which of the 
two processes was responsible, or to which degree. For the 
interpretation of radiocarbon dates of sites with an attri- 
bution to farming subsistence, one cannot find out whether 
demic or cultural diffusion was responsible for the appar- 
ent distribution in space and time of these sites. Or, put 
differently, the question on demic or cultural vectors may 
be not the most critical one for understanding the Neo- 
lithisation of Europe as a whole. 

What was the contribution of local adoption and inven- 
tion? From Figure [6] it is evident that for most regions 
conversion of resident foragers to farmers played a larger 
role than immigration. We quantitatively examined the 
relative importance of different sources (demic diffusion, 
cultural diffusion and local adoption or invention) to local 
Neolithisation in different model configurations (Tabic [TJ. 
Even in a scenario where demic diffusion is the only active 
process, migration as a source does not explain 100% of 
the agropastoralists in any focus region along the transect 
A-H but at most 85% (in region H), less in regions closer 
to the Mediterranean coast (e.g. 54% in region F). The 
local source (adoption and innovation) for an exchanged 
commodity like technology is in all configurations and fo- 
cus regions more important (70-90%) than migration or 
exchange. 



4-3. Independent agropastoralism 

Was regional exchange (via migration or trade) neces- 
sary at all for the onset of agropastoralism everywhere in 
western Europe? In Figure [3] (lower right panel) we show 
the result of a simulation where both exchange processes 
were suppressed, thus endogenous transitions to agropas- 
toralism become visible. The timing of the onset of agro- 
pastoralism around the Mediterranean Sea exhibits — next 
to the Levantine, Greek, and Moroccan centres which also 
appear in the reference simulation — many centres of hypo- 
thetical independent agropastoralism. This independent 
agropastoralism is solely predicted on the basis of suit- 
able environmental conditions (open vegetation type, not 
too cold) and internal development of sociocultural traits 
and demography; it corresponds to indigenist scenarios 
proposed in the older literature, for example for south- 



ern France (Geddes 1980 Courtin and Erroux 1974) or 



Greece (Theocharis and Bokonyi 1973 Winiger 1998). 



The indigenous agropastoralism hypothesis is, however, 
currently disregarded in archaeology because the genetic 
evidence points to the Near East as the centre for all 
Neolithic cultigens and nearly all domestic animals. For 
northern Europe, GLUES does not simulate the emergence 
of agropastoralism without the contribution of migration; 
these regions critically depended on the introduction of 
technologies and economies through the actual movement 
of people, commodities, and information. 

4-4- Model comparison and outlook 

In addition to the prior approaches to simulating the 
European Neolithic pavison et all 120061 120071 120091 |Ack 



land et al. 2007) which use geographic and topographic 



constraints for describing environmental heterogeneity, our 
model considers vegetation. Vegetation production is dir- 
ectly coupled to the carrying capacity and it determines 
the economic potential of a given environment. Already 
many regions, mainly around the Mediterranean Sea, have 
a high propensity for developing independent agricultures 
based on the palaeoecological background (Figure [3j lower 
right panel). We couple the diffusion rates of traits and 
migration of people not only to the background geography 
but also to the (evolving) technology. With these as- 
sumptions we can realistically reproduce the spatiotem- 
poral pattern of Neolithisation in greater detail than was 
done for the front speed of Neolithisation by |Davison et al.| 
( 20061 [20071) andlAckland et all (12007). 



Ackland et al. (2007) and more recently Patterson 



et al. (2010 for the Indian transition to agropastoralism) 



use the concept of converts to describe resident foragers 
which have converted to agropastoralism. We have shown 
that these converts may have played a larger role in the 
European Neolithisation than immigrant farmers. Our 
model provides additional insight into the processes re- 
sponsible for local adoption — often, a small share of in- 
troduced technology is sufficient to spark local invention 
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and trigger the transition. Alternatively, a few immig- 
rant farmers and the technologies and economic possibil- 
ities they carry along may suffice to stimulate the local 
transition. 

Our regional prediction for western Eurasia emerges in 
the context of a global simulation: not only is the subcon- 
tinental prediction embedded in the larger spatial scale, 
but every local transition occurs within the temporal con- 
text of preceding predominant foraging subsistence with 
continuous innovation and succeeding intensification peri- 
ods. While we have shown that the model realistically 
reproduced the European Neolithisation, where archeolo- 
gical data is plenty and most reliable, the model's spati- 
otcmporal consistency gives us confidence to draw conclu- 
sions about regions outside Europe in further studies. 

With the expected availability of more reliable palaeo- 
climate and palaeovegetation reconstructions from both 



models and data (Kutzbach et al. 2010 Gaillard et al 



2010), we expect to refine the large-scale biogeographic 



context of cultural evolution and the impact of local en- 
vironmental disturbances (Wirtz et al. 2010). We will 



then be able to assess better the degree to which the en- 
vironment determined the potential transition to farm- 
ing. This potential should, however, be interpreted in 



G. Ackland et al. (2007)'s way as providing a "historical 



null hypothesis. Its predictions can be taken as requir- 
ing no special explanation, and its failures can be taken 
as evidence of rare events that had significant and long- 
lived consequences". Numerical modelling of culture as 
a (natural) ecosystem may help to isolate the significant 
and non-deterministic events and concentrate our histor- 
ical interpretation on those events where culture was most 
emancipated from the environment. 

5. Conclusion 

We presented a spatially explicit mathematical model 
of the Neolithisation of western Eurasia from 8000 BC to 
3500 BC. Our model incorporates endogenous sociotech- 
nological dynamics, where culture is represented by the 
adaptation of characteristic population traits (technology, 
fraction of farmers, and economic diversity) and their in- 
teraction with demographics. The study resolved the spa- 
tial expansion of Neolithic culture via indigenous develop- 
ment, migration and information exchange and reproduced 
the chronology of agropastoral onset observed in field data 
across western Eurasia, particularly reproducing and ex- 
plaining the discontinuous speed of the 'wave of advance'. 

Our results encourage us to rethink possible indigenous 
centres along the northeastern shore of the Mediterranean: 
these might have not been able to develop since they were 
overrun by Near Eastern populations. Alternatively, the 
evidence for independent agropastoralism may have gotten 
lost in the admixture with the Fertile Crescent Neolithic 
package. According to our simulations, a north African 
contribution to the European Neolithic should equally not 
be discounted. 



The assessment of the relative importance cultural dif- 
fusion and demic diffusion in the model shows that either 
of these processes can explain the spatiotemporal pattern 
of agropastoral onset in Europe equally well. The phen- 
ology of the spatiotemporal pattern of agropastoral onset 
cannot discriminate the underlying process. Furthermore, 
even if only migration was considered and the diffusion 
of traits occurred only via immigrants, the prevalence of 
immigrant farmers in any of the emerging agropastoral re- 
gions was much less than the prevalence of foragers who 
adopted the agropastoral life style. To the long-lasting dis- 
pute between cultural and demic diffusionists our novel in- 
terpretation offers a balanced explanation of predominant 
adoption despite migration. Whether the adopting pop- 
ulation, however, did not also ultimately originate from 
Anatolia needs to be investigated by further archaeogen- 
etic studies. 
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